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Abstract 

A certain class of partial differential equations possesses singular solutions having discon- 
tinuous first derivatives ( "peakons" ) . The time evolution of peaks of such solutions is governed 
by a finite dimensional completely integrable system. Explicit solutions of this system are 
constructed by using algebraic-geometric method which casts it as a flow on an appropriate 
Riemann surface and reduces it to a classical Jacobi inversion problem. The algebraic structure 
of the finite dimensional flow is also examined in the context of the Toda flow hierarchy. Gener- 
alized peakon systems are obtained for any simple Lie algebra and their complete integrability 
is demonstrated. 

1 Introduction 

Camassa and Holm [1993] introduced and studied classes of soliton-type weak solutions, for an 
integrable nonlinear equation derived in the context of a shallow water model. In particular, they 
described the soliton dynamics in terms of a system of Hamiltonian equations for the locations of 
the "peaks" of the solution, the points at which its spatial derivative changes sign. In other words, 
each peaked solution, or peakon, can be associated with a mechanical system of moving particles. 
New systems of this type were obtained in Calogero [1995] and Calogero and Francoise [1996]. The 
r-matrix approach was applied to the Lax pair formulation of an n-peakon system by Ragnisco and 
Bruschi [1996], who also pointed out the connection of this system with the classical Toda lattice. 
A discrete version of the Adler-Kostant-Symes factorization method was used by Suris [1996] to 
study a discretization of the peakons lattice, realized as a discrete integrable system on a certain 
Poisson submanifold of gl(n) equipped with r-matrix Poisson bracket. 

In Alber et al. [1994, 1997, 1999] the existence of peakons was linked to the presence of poles 
in the energy dependent Schrodinger operators associated with integrable evolution equations. 
Namely, it was shown that the presence of a pole in the potential is essential for a special limiting 
procedure which allows for the formation of "billiard' weak solutions. By using the algebraic- 
geometric method, these billiard solutions are related to finite dimensional integrable dynamical 
systems with reflections. In this way, profiles of billiard weak solutions are associated with billiard 
motion inside quadrics with and without the presence of a Hooke's potential. The points of impact 
on the quadrics correspond to the peaks of the profiles of weak solutions of the nonlinear PDE. 
Billiard solutions include new quasi-periodic and soliton- like solutions, as well as peaked solitons 
with compact support. This method can be used for a number of equations including the shallow 
water equation, the Dym type equation, as well as A^-component systems with poles, together with 
all the equations in their respective hierarchies. 

In Section 2 of this paper we derive equations for the motion of the peaks of billiard solutions 
in the context of the algebraic-geometric approach. We construct solutions of these equations by 
connecting them to a flow on an appropriate Riemann surface which leads to a classical Jacobi 
inversion problem. 

In Section 3 we show that orbits of the Toda flow in sl(n) into which the peakons lattice is 
mapped by means of the Lax matrix introduced in Ragnisco and Bruschi [1996] admits a natural 
generalization to the case of an arbitrary Lie algebra g. Restrictions of the Toda flow in g to these 
orbits can be viewed as generalized peakons lattices. Using the methods developed in Gekhtman 
and Shapiro [1999] we give an intrinsic description of these orbits and prove complete integrability 
of the corresponding flows. Unlike in the sl(n) case, Chevalley invariants of g are not sufficient to 
ensure integrability in the general situation, and so we produce the requisite number of additional 
first integrals. Lastly, we explain how to construct Darboux coordinates for generalized peakon 
orbits and consider sl(n) and G2 as examples. 
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2 The Dynamics of the Peak Points 



As shown in Alber et al. [1994, 1999], quasi-periodic solutions of the shallow water equation 

U t + 3UU x = U xxt + 2U x U xx + UU xxx , (2.1) 
on the infinite line can be represented in the form 

U(x,t) = m + ---+iJt n -m (2.2) 
where //'s are solutions of the following systems of equations in x and t, 

Mi = "5— = Sign(^,m 2 i-i,m 2 i,Si) — ^j^— (2.3) 

and 



with 



Mi = = Sign(//j, m 2 i-i,m 2 i,di ) — v -, 

i = l,...,n, E = //H hMn, (2-4) 



2n+l 



r=l 



Here the constants m r depend on the initial condition for the equation ( |2.1| ), m = m r , and 

Sign(y, /, r, s) is a function on the Riemann surface, T = {w 2 = , which switches from —1 to 

1 and back each time y reaches I or r, endpoints of a cut on the Riemann surface. These ordinary 
differential equations (ODE's) for the fi- variables provide half of the equations of a finite dimensional 
Hamiltonian system. Notice that periodic solutions form a subclass of the quasi-periodic solutions. 
Also, soliton solutions can be obtained from quasi-periodic solutions by shrinking pair wize the 
roots m r of the spectral polynomial R{E). 

Peaked quasi-periodic solutions of equation ( |2.1[ ) can be constructed from solutions of (2.3) and 



(2.4) by using the limiting process mi — > (see Alber et al. [1994,1999]) and the trace formula 
(2.2). In what follows ODE's governing the time evolution of the peak locations are obtained and 
a connection with the canonical variables of Camassa and Holm [1993] is described. 

The x and t evolution of the fi variables can also be connected to a constrained motion of 
a particle on an n-dimensional hypersurface imbedded in M. n+1 which is parameterized by the 
constants m r , under the action of a harmonic force field (see Alber et al. [1999]). The limit 
mi — ► corresponds to "flattening" this surface along one direction which results in harmonically 
forced motion in a region of W 1 confined by an (n — l)-dimensional boundary. The particle collides 
with the boundary and moves along segments of an n-dimensional billiard. The reflection on the 
boundary causes, after using the trace formula (|2.2| ), the appearance of peaks in the corresponding 
PDE solutions of equation ( p.l|). 

On the level of (|2.3[) and (l2.4j) , the boundary of the billiard can be shown to be described by 
the following condition: fii(x,t) = 0. Reflection on the boundary is described by the sign switch 
from 1 to -1 (or from -1 to 1) of first derivative of fj,i in (2.3) and (2.4). 
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The Periodic Case. The main steps in deriving equations for the peak locations can be demon- 
strated by considering the following basic example. Let us assume n = 2, and mi = 0. Equa- 
tion fj,i(x,t) = defines a function of time x = q(t). An equation for the time dependence of q(t) 
can be found by differentiating (fj,i(q(t),t) = 0) with respect to time. We have 

^(Ml(g(t),*)) =Q=^(q(t),t)q+^(q(t),t). (2.6) 
The commuting x- and t-flows of [i\ are described by 



— = Signal , (2.7) 

OX {l\ - ^2 



dpi G - , J-^)\/C^i) 

— = Signal) * , 

at m - fi 2 

C 4 (/i) = (// - m 2 )(fJ. - m 3 )(fj, - m 4 )(fj, - m 5 ). 

Notice that the spatial derivative of fi± is not defined at x = q(t), so that in the above formula 
dfii/dx evaluated at (q(t),t) needs to be found in the limit as x — > q either from the left, or from 
the right. The sign uncertainty connected with this choice factors out from the final formula for q. 
This combined with ( |2.6D and the billiard condition at the boundary fj,i(q(t),t) = gives 

q = H 2 (q(t),t). (2.8) 

We also need an evolution equation for the quantity fj,2(q(t) , t) . Notice that from the trace formula 
( |2.2| ) and the expression for the boundary fii(q(t),t) = 0, this is the quantity that determines the 
amplitude of the peak U(q(t),t). The time derivative of fA 2 (q(t),t) is 

=^(Qm)q+^(q(t),t), (2.9) 



which used with the commuting x- and i-ffows equations for \i<i yields, in analogy with system (J2.8| 



d J t {t)=-^y)^(q + B l{ y)) 

= -Sign( y )^5M<j. (2.10) 

y 

Here we have introduced the notation 

y{t) = &t(q(t),t) 

and used again the definition of q(t) so that 

BiMq(t)) = Bt(y) = -fn(q(t),t) = 0. 

Combining this equation for y with the evolution equation (^^) for q gives the following system 
for the variables y and q, 



y = - Sign(y) (y - m 2 )(y - m 3 )(y - m 4 )(y - m 5 ) 

q = y (2.11) 
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Notice that the evolution for y is decoupled from that of q. The inversion problem corresponding 
to the equation for y determines a periodic function of t with nonzero average. The evolution of q 
is then given by the combination of a linear growth in time with slope given by the average of y, 
plus periodic oscillations. We also remark that for weak solutions of equation fl2.ip having isolated 
discontinuities in the first derivative, the jump condition yields 



and using the trace formula (2.2) this can be seen to coincide with the second equation in sys- 
tem (fPH) . 



The Genus 2 Quasi-periodic Case. We describe a billiard genus 2 quasi-periodic solution. The 
spectral polynomial in this case is of the 5th order and will be denoted as C^(fi). The boundary is 
again described by ni(q(t),t) = 0. We shall write 

2/1 = M2 (q(t),t) and y 2 = fi 3 (q(t) , t) . 

From the weak solution approach it follows that q = U{qi) = y\ + y 2 . We differentiate y\ and y 2 
to obtain 



2/1 



2/2 



<9/i 2 . d[i 2 



U*3(9,*)] 



dx 



dx 



dt 
dt 



Sign(fi 2 )^/C 6 (fi 2 ] 
^2^2 ~ M3) 



Sign(/i 3 )VC6(/^ 



(<? - Ma) 



Signal )y / C 6 (j/i 
2/1 - 2/2 



Sign(2/ 2 ) ^6(2/2) 



(2.12) 



^3(^3-/^2) ' 2/2-2/1 

This is a well-defined system associated with a genus 2 Riemann surface. The corresponding 
problem of inversion involves only holomorphic differentials and therefore the quantities yi,y 2 and 
q(t) can be expressed in terms of standard ^-functions on the Riemann surface. 



Peakon solutions. We now specialize the above formalism to the limiting case of soliton solutions 
of equation fl2,l|) on the real line, in which each pair m 2 i,m 2 i + i is taken to limit to the ctj, i = 
1, . . . , n. For just one /x-variable, using the trace formula U = fj, — a at the billiard boundary results 
in U(q(t),t) = —a. The single peakon solution determined by this procedure is therefore 

|a;+at| 



U(x, t) = —ae 

which is a traveling wave soliton-type solution. Camassa and Holm [1993] found 

U(x,t) = p(t)e-\ x - q W , 



(2.13) 



(2.14) 



with an additional link between p(t) and q(t) (Hamiltonian structure) which in the 1-peakon case 
results in p(t) being a constant and q(t) being a linear function. 



2-peakon solution. In this case, one obtains the following system that describes 2-peakon pro- 
files, 



d/J-i c . , , (fix - ai)(/zi - a 2 ) 
- Sign(/ii)- 



dx 

9/J.2 
dx 



Sign(/i 2 



(Mi - M2) 

t Qfg - Qi)(M2 - Qg) 
(M2 - Mi) 



(2.15) 
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where mi and M2 are evaluated between 01 and and a 2 and 0, respectively. The corresponding 
time-flow is 



9^1 1 \d 1 -°i)(Mi -02) 

- Signal - 



at 



Sign(/i 2 )-Bi(^ 2 



(Mi - M2) 
O2 - oi)(a*2 - a 2 ) 



dt ' ' ~ ' - (n-2 - //j ) 
and the first order polynomial i?i for equation (^[l]) gives 

#i(Mi) = -M2 + ai + a 2 , #i(M2) = ~Mi + a i + a 2 • 
Define the boundary by introducing functions qi(t) and 52 (i) such that 

Ml (<7l (*),*) =0, ^2(92 (*),*) =0. 

Define functions 

2/1 = /i2(?l (*),*), 2/2 = Ml (?2 (*),*)■ 
Differentiating Mi(9i(i) 5 £) and fi 2 (q 2 (t),t) i n (|2 ■ 1 7| ) results in 
^Ml n (dqi \ dm (dqi 



dfi 2 







<9a*2 / <ig 2 



+ Si(/i 2 ) 



9// 2 / d<?2 



2/2 + ai + 02 



cit dx \ dt J dx \ dt 

which coincides with the jump conditions for weak solutions 

~f~ = ^ (^i) = (Ml + M2 - ai - a 2 ) | ?i = M2(<?i) - ai - a 2 = yi - 01 - 02 
-^jr = ^(92) = (Ml + M2 - ai - a 2 ) | g2 = Mi (92) - «i - «2 = 2/2 - ai - a 2 



Finally, differentiate y\ and y 2 to find 



dyi _ <9/i2 rfQi ^M2 _ <9m2 ( dqi 



dt dx dt dt dx \ dt 
dy 2 d\x\ dq 2 d\i\ d\i\ ( dq 2 



dt 



+ 



dx dt dt dx \ dt 



— +^l(M2. 

+^i(mi: 



and so 



dyi 
dt 

dy 2 
dt 



Sign(yi)(yi - ai)(yi - a 2 ) 
Sign(y 2 )(2/2 - a 1 )(y 2 - a 2 ) . 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 
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Thus, the equations of evolution for y\ and y 2 decouple from those of q\ and q 2 . The decoupled 
equations ( |2.2lD can be solved first and q%, q2 subsequently determined from ( |2,20| ) by quadratures. 

It is interesting to examine the connection of the set of variables qi, yi, i = 1, 2, with the qi,Pi, 
i = 1,2, introduced by Camassa and Holm [1993]. By definition, the g's are the same in both sets 
(the (q,y) system does have a canonical Hamiltonian form). As to the y's, notice that 



U(qi) 
U(q 2 ) 



yi-ai-a 2 =Pi+ p 2 e' 



y 2 ~ai- a 2 =P2+ Pie 



(2.22) 



which provides a definition of the transformation of the y variables in terms of p's and g's. This 
together with ( [2.20 ) yields the first set of equations for the Hamiltonian system derived by Camassa 
and Holm [1993], 



at 

d -^= P2 +Pie-\^\ 
at 



(2.23) 



The constants of motion a\ and a 2 can be expressed in terms of p : s and c/'s via the first integrals 

P12 =Pi+p 2 = ~{a\ + a 2 ), 

and 



H 



12 



\{p\+vl)+PiP2e-^-^ = \(al + al), 



which are the total momentum and Hamiltonian for the (g,p)-flow, respectively. Using these 



expressions, the constants of motion a\ and a 2 can be eliminated from ( p.22|) and the explicit 
variable transformation of variables for yi, y 2 in terms of the (q,p) system is 



2/1 = P2 (e-^-^ 



-191-92 1 



y 2 =Pi(e 



1) 
1) 



(2.24) 



Differentiating these equations with respect to time, using d2T20| ), and the transforma- 

tion ( ^24]) itself, results in a system of equations for p\ and p 2 , which can be solved to yield 



= sgn(<? 2 - qi)pip 2 e 191 y2 ' 

at 

^ = sgn(<?i - q^p^e-^ 1 '^ 
at 



Equations ( 2.23| ) and ( |2.24| ) are precisely those that follow from H12 with a canonical Hamiltonian 
structure in the (p, q) variables. 

The case of three or more derivative-shock singularities Xi,yi, i = 1,2,... ,N > 3 proceeds in 
complete analogy with the case N = 2 above. Once again, the (/-flows decouple from those of 
the y's, however the equations in the system governing the y-flow are now coupled and it is not 
immediately obvious that this system is integrable. A closer inspection however reveals that the 
y-flow shares the same structure as that of the //-variables flow and is therefore integrable by a 
similar argument. 
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3 Generalized Peakons Lattices 



In this section we discuss the connection between the peakons lattice and the Toda flows. The fact 
that the peakons lattice can be realized as a special case of the Toda flows in gl(n) was discovered 
in Ragnisco and Bruschi [1996] and later used by Suris [1996] to introduce a discrete time peakons 
lattice. Let us review these results, concentrating, for simplicity, on a particular case of the peakons 
lattice with a Hamiltonian 

1 n 

h (p, <?) = 2 + ^2 PiP i ex P(^' ~ 1i) , (3-1) 

i=l l<i<j<n 

that corresponds to the case when particles qi are ordered as follows: q\ > q2 > ■ ■ ■ > q n - It turns 
out that (3.1) possesses a Lax representation with the Lax matrix 



L = L(p,q) = ^T PiEu + ^2 VPiPj ex P ( 2 (H ~ *)) ( E v + E ii)i ( 3 - 2 ) 

i=l l<i<j<n ^ ' 

and an auxiliary matrix equal to the half of the skew-symmetric part of L, the peakons lattice thus 
being a restriction of the Toda flow to the set of matrices of the form (|3.2| ). It was shown that 
this set forms a Poisson submanifold w.r.t. the associated r-matrix bracket on gl[n). Complete 



integrability of (3.1) is provided by the Chevalley invariants of gl(n), i.e. spectral invariants of L. 

In what follows we give a uniform Lie-algebraic procedure which allows us to construct peakon- 
type orbits in any simple Lie algebra. On each of these orbits the Toda flow is shown to be a 
completely integrable system. Notice that for proving this we provide additional first integrals 
which supplement the Chevalley invariants of the algebra. We also indicate a convenient way of 
constructing Darboux coordinates for orbits under consideration. 

Let us first recall the Hamiltonian formalism for the generalized (symmetric) Toda flows (cf. 
Kostant [1979], Goodman and Wallach [1985], Reyman and Semenov-Tian-Shansky [1994]). Let g 
be the normal real form of a simple Lie algebra of rank r, G a corresponding Lie group and f) a 
Cartan subalgebra of q. We denote by $ the root system of q, and by <£ + (resp. <3? _ ) the set of all 
positive (resp. negative) roots. We also fix a Chevalley basis {e a , a € hi, i = 1, . . . , r} in Q . 
All properties of the root systems and Chevalley bases that we need can be found in Humphreys 
[1980] and Onishchik and Vinberg [1990]. 

Consider a direct sum decomposition q = t + b + , where b + is the upper Borel subalgebra and 6 
is the maximal compact subalgebra of q. The dual space b?j_ of b + can be identified with the space 
S of "symmetric" elements of q that have a form 

r 

L=} J bjhj + a a (e a + e_ a ) . (3.3) 

j=l ag<I> 

S is an annihilator of I w.r.t the Killing form ( , ). 

The pull-back of the Lie-Poisson bracket on b* + equips S with the Poisson bracket 

{fuh}s{L) = i(L,[7r + V/i(L),7r + V/ 2 (L)]) , (3.4) 

where gradients are defined w.r.t. the Killing form and 7r+ is a projection on b+ parallel to 6. 
Equations of motion of the generalized Toda flow on S are generated by a Hamiltonian H{L) = 
\{L, L) and have the Lax form 

L=[L,±ir + (L)] = [±w-(L),L] (3.5) 
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where 7r_ = Id — tt + . The bracket ( |3.4| ) is a restriction to S of the so called r-matrix bracket on g: 

{fx, h}r{X) = ~(X, [7T + Vh(X), 7T + V/ 2 (X)] - [7T_V/l(X), 7T_V/ 2 (X)]) . (3.6) 

Symplectic leaves of the Poisson manifold (S, {,}s) coincide with the orbits of the coadjoint 
action of the upper Borel subgroup B + of G: 

Ad* b (L) =7r 5 Ad 6 -i( J L), (3.7) 

where its is a projection on S parallel to b+. On every indecomposable symplectic leaf, restrictions 
of the Chevalley invariants I\ = H, I2 , ■ ■ ■ ,I r of q form a family of independent Poisson commuting 
integrals for the Toda flow. This family is maximal, however, only for a few distinguished orbits, 
that were classified in Goodman and Wallach [1984] and Perelomov and Kamalin [1985] 

To describe a particular type of coadjoint orbits, which in the sl(n) case will be shown to be 
associated with the peakons lattice, we need the following definitions introduced in Gekhtman and 
Shapiro [1999]. Let m be the maximal positive root and 

Define 

0' = Ker ad Cm nKer ad e _ m 
F = Span{e a : a G and (m, a) / 0} (3.9) 
F = Span{F, h m } 
V = F-R(e. m ) . 

Then q' is a semisimple subalgebra of q and F is a Heisenberg subalgebra of g, i. e. V is spanned 
by root vectors e ai , e- m - ai , i = 1, . . . N , such that 

[^Oij^aj] — m— otii m— Oj] — [^— mi &an\ — [fi-mi 6-m-a,] — 0, 

\&a.i i 6— m— Oj] — Ci5^€— m , , (3.10) 

where Cj are some positive constants. Note also that 

[h m ,e ai ] — e ai , [h m , e_ m _ Q!i ] — e_ m _ Q , i , [/im)fi±m] — i2e-i- m , (3.11) 
Therefore ([O]) and ( |3.11 ) yield the following Poisson brackets 



{yo,x }s =x , {yo,Xi}s = 2 Xi > {yo,yi}5 = \v^ {Vi, x i}s = 7^0, i = l,...N, (3.12) 
for the linear functions 

Xi = (L,e- ai ), yi = —(L,e m+ai ), x = (L,e m ), y = (L,h m ) . (3.13) 



All other brackets between functions ( B.13 ) are zero. 



Let ho € f) be orthogonal to h m w.r.t the Killing form and let Om,ho be the coadjoint orbit of 
B + through (e_ m + h + e m ). 
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Theorem 3.1 Q m ,h can be parameterized by elements of F : 







m,ho 



L = TT S [( 



1 



MO 

where element Q G F- is decomposed as 

C = x(C)e_ 



ad em v((),v(()]) :(GF.,x(C)>0 



(3.14) 



y{()h m + v(() 



and v(Q G V. 

The Toda flow ( \3. 9J is completely integrable on O mt h with a maximal family of functions in 
involution given by the coefficients of polynomials 



(3.15) 



where F are the Chevalley invariants of i 



The proof makes use of the so-called 1-chop map, a Poisson map from q to g', which was introduced 
for sl(n) in Singer [1990, 1993] and generalized for an arbitrary Lie algebra in Gekhtman and Shapiro 
[1999]. This map and the functions ( |3.15 ) play a key role in the proof of the complete integrability 
of the Toda flows on generic orbits (see Deift et al. [1986] and Ercolani et al. [1993] for the sl(n) 
case and Gekhtman and Shapiro [1999] for a general case). 

It also follows from Theorem 2.1 that for constructing Darboux coordinates on D m ^ , it is 
sufficient to construct Darboux coordinates for a rather simple Poisson algebra gener ated by (gig ). 

The following example demonstrates the connection with the peakons lattice. 

= E n i,h m = En — E nn , 



JJ 



m — 


Ein, e 














T 


-v . 


4 



: x,y &R,V\,V2 & 



D n-2 



Since Tr(L) is invariant under the action (3.7), we can define S to be the space of all symmetric 
matrices, not necessarily with zero trace. One can choose ho to be any diagonal matrix of the form 
ho =diag(«;, d, k), where d is a (n — 2) x (n — 2) diagonal matrix. It follows from ( 3TJ ) that 



D 



m,ho 



L 



y + « 



X 



X 

K-y 



: x > 0, y G R,vx,V2 G 



on— 2 



(3.16) 



If d = 0, an open set { — k < y < k\ vu > 0, V2% > 0, i — 1, ... ,n — 2} C £)m,ho admits Darboux 
coordinates pi and qi defined as follows: 



Pi = y + k, p n = K-y, p i+ i 



V\iV2i 



x 



1,... ,n-2, 



and 



exp(q n 



K 2 _ y2 



, exp(qi - qi) 



XVu 



(y + K)v 2 i 



2,... ,n-l. 
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With this parameterization matrix L in ( |3.16| ) coincides with the Lax matrix fl3.2j). Note that 
Pi + • • • + p n = const on £> m ^h an d that coordinates qi are defined only up to a translation. Note 
that a different choice of Darboux coordinates for D mj o was proposed in Kamalin and Perelomov 
[1985] as an example of the general construction of canonical coordinates on coadjoint orbits. 

Example 2.3. Let g be an exceptional algebra of type G 2 - Denote the short and long simple 
roots of Q by v\, v 2 resp., and root vector corresponding to a negative root a = — (iv\ +3^2) by dj. 
Then F is spanned by vectors e_ m = 632 and en, i = 1, ...4. Thus, for any ho the orbit O m ,ho 
is 6-dimensional, while there are only two independent Chevalley invariants. This illustrates the 
necessity of using functions ( |3.15 ) to establish complete integr ability. 

Put ho = and choose Darboux coordinates qi,Pi for Poisson algebra ( |3.12| ) with N = 3 in the 
form: 

yo xi 6yi y 2 2x 2 
qi = x , pi = , 92 = -^-f=, P2 = —p=, q3 = f=> P3 ~ 



xo Z^/xo v^o Vxo y/Xo~ 

Then Theorem 2.1 gives the following integrable polynomial Hamiltonian quadratic in momenta: 

H(p, q) = ^( P3 q2 + e^ 2 ) 2 + ^(p 3 q 3 + p 2 q2? + ^(pI + 12 pI + 3p?) + ?? + Mhl + 3^ (3.17) 

A detailed analysis of the system generated by ( p. 17 ), as well as integrable systems associated 
with peakon-type orbits in other simple Lie algebras will be given in a forthcoming publication, 
Alber and Gekhtman [1999]. 
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